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o 

Energy spectra and electromagnetic transitions of nuclei are strongly depending from 
the correlations of the bound nucleons. Two particle correlations are responsible for the 

■ scattering of model particles either to low momentum- or to high momentum-states. The 
q ' low momentum states form the model space while the high momentum states are used to 

Q calculate the G-matrix. The three and higher order particle correlations do not play a role 
in the latter calculation especially if the correlations induced by the scattering operator are 

■ of sufficient short range. They modify however, via the long tail of the nuclear potential, the 
' Slater determinant of the A particles by generating excited Slater's determinants. In this 

work the influence of the correlations on the level structure and ground state distributions 
. of even open shell nuclei is analyzed via the boson dynamic correlation model BDCM. The 

' model is based on the unitary operator e s (S is the correlation operator) formalism which 

■^j- | in this paper is presented within a non perturbative approximation. The low lying spectrum 

QO , calculated for 6 Li reproduce very well the experimental spectrum while for 6 He a charge 

radius slightly larger than that obtained within the isotopic-shift (IS) theory has been cal- 
£SJ ' culated. Good agreement between theoretical and experimental results has been obtained 

t-H , without the introduction of a genuine three body force. 

\Q '■ PACS: 21. 10. Ma, 21.60.-n, 21. 10. Ft, 21.60.Fw 

O ; 

T ^ \ §1. Introduction 

Correlation effects in nuclei have been first introduced in nuclei by Villars, 1 ^ who 
proposed the unitary-model operator (UMO) to construct effective operators. The 
method was implemented by Shakin 2 ) for the calculation of the G-matrix from hard- 
core interactions. Non perturbative approximations of the UMO have been recently 
applied to even nuclei in Ref. 3) which here is treated in more detail. The basics 
formulas of the Boson Dynamic Correlation Model (BDCM) presented in the above 
quoted paper have been obtained by solving the n-body problem in the following 
approximations : 

a) The n-body correlation operator is separated in short- and long-range compo- 
nents. The short-range component is considered up to the two body correlation 
while for the long range component the three and four body correlation operators 
have been studied. The extension of the correlation operator to high order diagrams 
is especially important in the description of exotic nuclei (open shell). In the short 
range approximation the model space of two interacting particles is separated in two 
subspaces: one which includes the shell model states and the other (high momentum) 
which is used to compute the G-matrix of the model. The long range component of 
the correlation operator has the effect of generating a new correlated model space 
(effective space) which departs from the originally adopted one (shell model). The 
amplitudes of the model wave functions are calculated in terms of non linear equa- 
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tion of motions (EoM), 

b) the n-body matrix elements are calculated exactly via the Cluster Factorization 
Theory (CFT), 

c) by linearizing the systems of commutator equations, which characterize the EoM. 
The generalized linearization approximations (GLA) includes in the calculation pre- 
sented in the paper up to the (3plh) effective diagrams. The linearized terms are not 
discarded but provide, as explained later in the text, the additional matrix elements 
that convert the perturbative UMO expansion in an eigenvalue equation. 

Within the present treatment of the correlation operator one generates in the 
n-body theory not only the ladder diagrams of Ref. 4) but also the folded diagrams 
of Kuo . 5 ) 

In this paper the BDCM model is applied to calculate the influence of the corre- 
lations on the energy spectrum of 6 Li and on the charge distributions of 6 He and 6 Li. 
The motivation of these calculations relays from one side in the study the effect of 
the correlation operator on the theoretical charge radius of 6 He and from the other 
side in investigating the variation of the charge distribution and magnetic moment 
of the ground state of 6 Li under the variation of the adopted model space. 

The value obtained for the charge radius of the correlated 6 He is slightly bigger 
than the radius calculated in other theories 6 )' 9 ) and that derived within the isotopic- 
shift IS theory. 10 ^ A charge radius which agrees with the radii calculated in the 
Refs. 6, 9) and those calculated in the cluster models of Refs. 11-13) is on the other 
hand obtained by considering only two protons in the Is i . This non correlated radius 

2 

agrees also with the radius derived at Argonne within the IS theory. 10 -* Correlations 
have therefore the property to increase the charge radius of 6 He as observed for the 
isotopes of Lithium. 

The calculations performed in Ref. 14) for the charge radii of the lithium iso- 
topes, although in good agreement with those measured at GSI-TRIUMF 15 ) and 
analyzed with the help of Ref. 16), are always slightly larger than those measured. 
For the stable isotope 6 Li the calculated radius agrees with the value obtained with 
the electron scattering experiments of Ref. 17). However, the charge radii calculated 
in the IS theory could also depend on the nuclear correlations. The consideration 
of the microscopic correlations as presented in this paper will generate a new eval- 
uation procedure for the Mass Shift (MS) and the Field Shift (FS). As result, both 
quantities could be evaluated within non-perturbative methods which include the 
nuclear effects. 

The importance of the correlations in the evaluation of the FS has been already 
pointed out in Ref. 19) where we calculate the FS of 7 Li and show that the departure 
from a point nuclear approximation is rather a big effect. Additionally the higher 
order cross term contributions of Ref. 18) need to be considered. A direct comparison 
between the calculated and the measured charge radii should be therefore performed 
after an accurate analysis of these two correcting factors. 

Theoretically the effect of the correlation on the distributions of medium-heavy 
nuclei has been already performed in Refs 20-22) within a phenomenological cor- 
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relation model. In light nuclei, however the calculation of distributions performed 
by the theoretical models of Refs. 6, 7, 8, 9) has been done in terms of non corre- 
lated particles. In Refs. 6, 7) and references therein quoted large-basis no-core shell 
model have been used. In this model the e s method is considered up to the "two 
body cluster approximation". Within this approximation the effective interaction 
obtained contains no hole state. This will be however not the case by expanding the 
e s method to higher cluster approximations as done in the present paper. Unstable 
nuclei have been described by the antisymmetrized molecular dynamics (AMD) of 
Ref. 8) neglecting correlations. Quantum Monte Carlo calculations have been per- 
formed by using realistic nuclear Hamiltonians that fit nucleon-nucleon scattering 
data in light nuclei. 9 ) Good results have been obtained for structure calculations, 
but the model can not introduce correlated wave functions in the calculation of the 
distributions. 

Another important argument for the consideration of correlation effects comes 
from the analysis of the magnetic moments. The magnetic moment of the ground 
state of 6 Li calculated in a model base in which the hole is confined to the lsi is 

2 

smaller than the experimental value. Only with the use of a large configuration base, 
which includes the spin-flip component lp^lpi proposed by Arima 23 ) we obtain a 

2 2 

very good value for the magnetic moment. The large configuration base has also the 
effect of decreasing the energy of the second 1~ level. 

§2. Theory of Correlated Two Particle Systems 

The effect of the correlations between nucleons in open shell nuclei is investigated 
within a system of coupled commutator equations, which via the GLA Ansatz, is 
converted in an eigenvalue equation. These describe a situation in which valence 
particles and core-excited states are coexisting. The advantage of this model is to 
provide larger effective configuration spaces (see Appendix A) than that used in the 
shell model calculation, to include exactly the Pauli principle in the coupled spaces, 
and to generate correlated solutions for the n = paired systems. 

We start by computing the commutator of the valence particles with the nuclear 
Hamilton's operator. In this calculation we retain the linear (shell model) terms 
with dimension n and non-linear (valence pairs coupled to one particle-hole pair) 
excitations with dimension n+V. The 1' denotes one particle-hole pair. In the next 
equation we have then to take care of the non linear terms derived in the first step. 
The commutator of the Hamilton's operator with the n + V excited states generates 
the coupling of the n-valence pairs to two particles-holes (2p-2h) excitations with 
dimension n + 2' . The successive equations are then characterizing the commutators 
which involve valence particles coupled to an increasing number of particle-hole pairs 
n + n' . It is worthwhile to remark that the obtained system of commutator equations 
is similar to the chain of equations one derives within the Green's function dynamics 
of Ref. 9). The introduced computational steps describe the mixing of the shell model 
states to core excitations with an increasing degree of complexity which will find 
applications in the calculation of the structures of exotic nuclei. This commutator 
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chain is suitable to be solved perturbatively by inserting the ra-th commutator in 
the (n — l)-th commutator, (n — l)-th commutator in the (n — 2)-th commutator , 
. . . the second commutator in the first. Within this perturbative approach one defines 
effective Hamiltonians of the model which, due to the increasing degree of complexity, 
are not easily solvable. Much simpler solutions to the commutator equations may, 
however, be obtained in the BDCM model. We start by remarking that in the study 
of low lying excitations of the n-body systems the higher order components of the 
wave functions, which involve n valence - and (2p-2h) core-excitations are poorly 
admixed in the model space and can be linearized. Within this approximation, the 
model commutator equations are suitable to be restricted to a finite space. The 
linearized system of the commutator equations is then solved exactly in terms of the 
CFT which calculates the n-body matrix elements in an expedite and exact way. 

In the following we illustrate the method by considering two valence particles. 
Following Ref. 2) we calculate the effective Hamiltonian by using only the S 2 corre- 
lation operator obtaining: 

H eff = e-^He iS * = £ a/3 Hi|/?)4 a,3 + E^^l^l^aU^ 
= Y. a p{a\t\j3)aia p + Ea^^M^) ^ ^ 

where v\ 2 is the long component of the two body interaction (note that the v[ 2 is in 
the following equations simply denoted as v). The is the two particle correlated 
wave function: 

9 a ft = e iS2 <P af 3 (2-2) 

In dealing with complex nuclei however the (Si, i = 3 ■ ■ ■ n) correlations should 
also be considered. 

The evaluation of these diagrams is, due to the exponentially increasing number 
of terms, difficult in a perturbation theory. 

We note however that one way to overcome this problem is to work with e t ( Sl+S2+S3 ~ { 
operator on the Slater's determinant by keeping the n-body Hamiltonian unvaried. 

After having performed the diagonalization of the n-body Hamilton's operator 
we can calculate the form of the effective Hamiltonian which, by now, includes cor- 
relation operators of complex order. 

We write the two particle states in second quantization by discarding for sim- 
plicity the isospin quantum numbers: 

$2 P — A\( ai J)\Q) = [a] ia ] 2 ] J M \0), (2-3) 

where the operators a] i a] 2 create two coupled particles in the open shells and we 
analyze the structure of the particle dynamics, generated by the correlation operator, 
via the following commutator: 

[H,A\( ai J)]\0) = [(^e Q aK + ^(a/?| U (r)|75)4aJa (5 a 7 ),(a] i a] 2 ) J ]|0). 

a 0/375 

(2-4) 

In order to have a compact index definition we have introduced: 

"i — ► jih (2-5) 
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By using some operator's algebra and by including in the results linear and nonlinear 
terms we calculate: 

[H,A\( ai J)}\0) = J2tt(2p\2p')4(PiJ)}\0)+ E n(2p\3plh)4({3 2 J[4J)}\0). 

(2-6) 

In Eq. $FQ$) the A\ J) operators are those of Eq. ([13]) and the Aj,(/3 2 J( J' 2 J) are 
defined below: 

$zpih — 4((3 2 J[J'J)\0) = ((a],a],)^(a],a,,)^)^|0). (2-7) 
In Eq. (|2-7p we have used the additional convention: 

ft — ibsisii (2-8) 

and we have associated: 



J[ to the coupling of j[j 2 
J' 2 to the coupling of j 3 j'4 



(2-9) 



Having extended the commutator as in Eq. (|2-6p . we have also to calculate the 
commutator equation for the A\{a 2 J\J 2 J) operators as given below: 

[#,4(a 2 JiJ 2 J)]|0> 

= E ^(3pl/i|3p'l/i')4(/ ? 2«/i'/2 J )l ) + E ^(3pl/i|4p2/i)^(/3 3 J(J 2 J^J)|0), 

(2-10) 

where we have introduced the (4p-2h) wave functions defined below: 

$4 P 2h — 4(ft «^)|0) = (((a] i a^^(a],a i /) J i) J "(aJ,a ii ) J s)J|0) J (2-11) 
and where we have consistently extended the definition given in (|2-5l2-9[) : 

ft — iii^Vsie (2-12) 

with: 

J[ associated to the coupling of j[j 2 

J 2 associated to the coupling of (2 -13) 

J3 associated to the coupling of 

In the definition of A\{j3^ J[ J 2 J3 J) the coupling of J{ to J' 2 to J12 has been dis- 
carded from the notation. In Eqs. (|2-6|2-1U|) the fi's are the matrix elements of the 
Hamilton's operator in the model wave functions. The next step would be then the 
computation of the commutator of the Hamiltonian with the A3 (ft J[ J 2 J 3 J) oper- 
ators. Here we linearize these contributions by considering that in the study of the 
low energy spectrum and in the calculation of ground-state correlated distributions 
the ^-3 (ft J'i J 2 J 3 J) terms are poorly contributing. The linearization is performed as 
in Ref. 25) by applying to the (4p2h) terms: 

E {aPHrWaypas^ifofMJ) (2-14) 
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the Wick's theorem and to discard the normal order terms. Within this linearization 
approximation we generate from the commutator equations of Eq. (|2 6I2-10[) non 
perturbative solutions of the EoM, i.e.: the eigenvalue equations for the mixed mode 
system: 

[H,A\( ai J)]\Q) = ^(2p|2p')4(/W)|0) + Y, n(2p\3p'lti)Ai(f3 2 J[j!,J)\0), 
Pi p 2 J[J 2 

(2-15) 

and 

[H,A\(u 2 J x J 2 J)}\Q))= Y n ( 3 P lh \ 2 P) A i(M\°) + Y V(3plh\3p'lti)A\(l3 2 J[J' 2 )\Q). 

Pi P2J[J' 2 

(2-16) 

Within the application of the GLA approximation we convert Eqs. (|2-6I2-I0p in an 
eigenvalue equation for the configuration mixing wave functions (CMWFs) of the 
model. In fact, the linearization provides the additional matrix elements necessary 
to write the following identity: 

Q{3plh\3p'lh') = Ud^ajMrMfM), (2-17) 

and to introduce the off-diagonal matrix elements which couple the (2p) to the (3plh) 
subspaces. Now, by writing Eqs. (|2-15|2-16l) in the following matrix form: 

/ [H,A{( ai J)}\0) \ = f E 2p + Q{2p\2p') f2(2p\3p'lti) \( 4(M|0) 

I [H,Al(a 2 JiJ 2 J)]\0) ) V n(3plh\2pf) E 3plh + f2(3plh\3p'lti) J I A\(p 2 J x J 2 J)|0) 



and by multiplying to the left with: 



((pi(aiJ) 
(0\A 2 (a 2 JiJ 2 J) 



(2-18) 
(2-19) 



we generate the eigenvalue equation for the dressed particles: 

' E 2v + {A x {a x J)\v(r)\A\(p x J)) {Ai(a\J)\v(r)\A\(f5 2 J' l J 2 J)) 



E 

PiPij[j' 2 



{A 2 {a 2 J x J 2 J)\v{r)\A\W)) E 3plh + (A 2 (a 2 JxJ 2 J)\v(r)\ A\ (/3 2 J[ J 2 J) ) 



XiWiJ) \ _ rp ( Xl(ai^) 



X2(M4J) J E \ x2(a 2 J x J 2 J) J |0) - 

(2-20) 

In Ea. (l2^20ll E 2p = ef/+ef 2 F and E 3plh = ef/ '+ef/ +eff -ef/ are the Hartre-Fock 
energies (see Appendix C) while the x' s are the projections of the model states: 

\$l P ) = Xi(«i^)Ai(aiJ)|0) + X2(a 2 J 1 J 2 J)A 2 (a 2 J 1 J 2 J)\0) (2-21) 

to the basic vectors 2p, 3plh. To conclude, although the (4p-2h) CMWFs are not 
active part of the model space, they are important for structure calculations. On 
may therefore associate the GLA approximation to a parameter which describes the 
degree of complexity of the model CMWFs (the method used to define the CMWFs 
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is given in the Appendices A, B, C). Within the first order linearization we obtain 
the EoM for the shell model while within the second and third order linearization 
approximations we derive the EoM of valence particles coexisting with the complex 
particle- hole structure of the excited states. 

In this paper we solve Eq. (|2-20j) self-consistently (see Appendix D). The solu- 
tions for the first iteration step are obtained by diagonalizing the eigenvalue equa- 
tion (|2-20p . The first step of the iterative method generates the dynamic amplitudes 
for the two dressed particles, i.e. two particles coexisting with the 3plh struc- 
tures. With the calculated eigenvectors we recompute then the matrix elements 
{hj2\v{r)\hfahh) and {hhhk^ir^j'^j'^ and we diagonalize again the eigen- 
value equation. The iterations are repeated until the stabilization of the energies has 
been reached. 

To be noted that since we are working in coupled particle-particle and particle- 
hole bases we need both particle-particle and particle-hole matrix elements (these 
terms where set to zero in the original work of Briickner 4 )) in orders to diagonal- 
ize the eigenvalue equation. The calculation of these matrix elements is, however, 
complicated by the number of terms which have to be evaluated in order to solve 
the introduced iterative equations. As in Ref. 3) one solves this problem by using 
the cluster factorization theory CFT which provides a quick and exact numerical 
method to perform the calculations of the matrix elements. The starting point of 
the CFT theory is to expand the CMWFs base introduced for (3plh) states in terms 
of cluster factorization coefficients (CFC) denoted in the following U and V: 

|3pl/i) J =4(a 2 JiJ 2 J)|0) 

+ 5 uf' 1 \a 2 JiJ2\}a j J r a J J s )\[B^a J J r )A\(a j J s )] J \0) 

JrJsOj (2-22) 

+ 5 uf' 1 \a 2 JiJ2\}a j J r a J J s )\[a j J r a j J s } J \0), 

where the (a 2 ) coordinates of the (3plh) model states have been expanded in terms 
of active (passive) particle-particle and passive (active) particle-hole coordinates as 
given below: 

"2 — ► £j(Jr)(active particle-particle) ej(J s ) (passive particle-hole) , . 

+aj(J r ) (active particle-hole)aj (J s ) (passive particle-particle) 

In Eq. (|2-22p the over-script on the left of the V and U indicates the total number of 
pairs (2(two pairs)+l=5) while the over-scripts on the right the number of particles 
and holes. There we have also introduced the creation operators of a particle-hole 
pair: 

(ph) B\( ai J)\0) = [ot lOia ]^|0) (2-24) 



8 



(M. Tomaselli et al.) 



and postulated that the sum over (i) is to be extended over the six combinations 
(partitions) of particle-particle pairs formed by the (3plh) model space and (j) over 
the three combinations of particle-hole pairs 

OL2 ► £j£j OijOLj . . 



Within this convention we reproduce exactly the the (3plh) CMWFs. In Eq. (|2-22p 
the sum over (oij) and (e^) has not been given explicitly because these indices are 
complementary to the (04) and (ej) in the sense of Eq. f)2-25|) . The sum over J r , J s 
take care of the fact that by calculating the (3plh) matrix elements in the traditional 
way, the coupling of the j's of the (2p) is not the same as in the coupling postulated 
for (3plh) wave functions. The definition of active and passive components separate 
the matrix elements of the interaction in particle-particle and particle-hole matrix 
elements. 

It has to be remarked that, due to the used iterative method, the final results 
are almost independent from the initial choice of the interaction because after the 
first iteration the two body potential is modified by the contributions of both types 
of matrix elements. 

Matrix elements calculated by using the first term of Eq. (|2-22p are related to the 
Bruckner theory, 4 ) while those evaluated with the second term have been considered 
in the folded diagram theory of Kuo. 5 ) 

In order to calculate the C¥C-^V j^ ,1 \a2J\J2\}^iJr^i J s) we introduce 3 ) the fol- 
lowing operator: 

<(2) = [(4a}) J «(o i o i 0^ (2-26) 

which destroys and creates a particle pair and we evaluate his effect on the 3plh 
wave-function. In Eq. (|2-22p the two in the parenthesis indicates that the operators 
are working on the space spanned by two boson (two pairs). We calculate: 



<(2)4(a 2 JaJ 2 JM)|0) = C^^Coei! • ([a\4] J '[a] 3 a jA ] J >) J M 
+Coef 2 • ([a\a) 2 ]^[a\,a H ]^y M ,+Coeh • {[a)^^^)^)- 



(2-27) 



The Coefx^Aj = 1,2,3 are given in the appendix E. By multiplying Eq. (|2-27|) to 
the left with (0^2 (02^1^2) we obtain the following equation: 

(0|A 2 (a 2 JiJ 2 JM)7r^ fc (2)4(q 2 JiJ 2 JM)|0) _ JkJ 

(Coefi+Coef 2 +Coef 3 ) ~~ U Mm k M' K > 

where we have introduced the Clebsch-Gordan coefficients Cfjf^ M , . Eq. (|2-28p de- 
fines the matrix elements of a unitary operator, i.e.: by applying the unitary operator 

7T fc (2) 

Coefi+Coef 2 +Coef3 on * ne °f a wave function of the type \3>plh) J M we reproduce 
the same wave function in a rotated frame. It follows that the structure defined 
by Eq. (|2-28|) is that of a full linear group in (2J+1) dimension and of its unitary 
subgroup t/2j+i- The calculation of the transformation coefficients is now reduced 
to the construction of the CoefA coefficients. In order to demonstrate this, let us 
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now operate with the Casimir's operator of the group (7r^ fe (2) <g> (^m k (^))^)o on the 
left side of Eq. (12^221) . We obtain: 

(yl 2 (a 2 J 1 J 2 J)||(<(2)®(vri fe (2))t)0||4(a 2 J 1 J 2 J)) = ^Coef Ai Coef Aj (2-29) 

where the double bar matrix elements have been introduced because the CFC are 
independent from the m's quantum numbers and where the subindex Aj classifies the 
three different partitions spanned by the particles in the (3plh) wave functions. On 
the other hand considering that: 

A; 

« fc (2) ® «(2))t)° = ® ((vr °(l^))t«(l,z))t)0 (2-30) 

i=l 

where the sum is running over all the possible partitions, we calculate by using 
Eq. (I2^22j) : 

(A 2 (a 2 J 1 J 2 J)\\(7T k mk (l,i)7r° (l,i)^ (( 7 r (l,z))t( 7 r^(l,i))t)0||4(a 2 J 1 J 2 J)) 
= 5 y) (3 ' 1) (a 2 J 1 J 2 |}A i e i J r e i J s ) 5 yj 3 ' 1) (a 2 JiJ 2 |}A i e i J r elJ s ). 

(2-31) 

In Eq. (|2-31|) we have introduced the unit operators 7Tq(1,«) defined by: 

( J B 1 (e J j;)|(7r (l,i))tvr (l^)| J B 1 t( e - jJs )) = S^Sj.j,. (2-32) 
By equating Eq. ([2 T 29l) and Eq. ([2 T 3T]) we obtain: 
5 F) (3 ' 1) (a2JiJ2|}A i eiJ r elJ s ) 5 yJ 3 ' 1) (a 2 JiJ2|}A i e i J r e i J s ) = ^ Coef Xi Coef Aj . (2-33) 

Eq. (|2-33p is given in matrix form in the Appendix E. Note that from the diagonal- 
ization of the Eq. (|E-4p we derive the square of the CFC. In order to calculate the 
CFC we introduce the same phase convention used to define the Clebsch-Gordan 
coefficients. Analogously also for the operator 

u k mk (2) = l(a\a J ) J <(a] / a l ,y>] k mk (2-34) 

which destroys and creates a particle-hole pair. The normalization factors for these 
operators (c^; i = 4,5, 6) are given in the appendix D. A numerical example for the 
calculation of these coefficients is given in Appendix F. By using Eq. (|E4l) and the 
evaluated CFC, we can now calculate the matrix elements of the nuclear interaction 
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in the (3plh) CMWFs. We derive: 

(3plh\v(r)\3p'lti) = (A 2 (a2JiJ2J)\v(r)\A^(P2J[J2J)) 

5 ^ J t (a2Ji«/ 2 |}A i e i J r e,J s ) 5 y 7 (/3 2 J(^|}A J e j j;e- i ^) 

Aj Xj tit j J r J s J !f J $ 

{[KeiJreiJsYH^WX^J'^J 1 ^) 

+ Yl 5 Ul(a2JlJ2\}^iaiJraiJ s ) 5 Uj((32J[J2\}H a j J rajJ' s ) 

([^iaiJ r aiJ s } J \v(r)\[^ja.jJ' r a.jJ' s ] J ) (2-35) 
Yl 5 ^j t (a 2 JiJ2|}A i e i J r e i J s ) 5 y J (/3 2 J{^|}A i e,j;e i J s ) 

Aj Xj Ei^-jJrJ^Js 

+ Yj 5 U\(a2JlJ2\}ViaiJrUiJs) 5 Uj(/32J[J2\}fJ'jUjJrajJ s ) 

Hi/J,jaiajJ r JrJa 

{^i\v(r)\fijaj} J a r . 

In Eq. (|2-35p we have two types of matrix elements, the particle-particle matrix 
elements which, by choosing one Aj and [i{ partitions, can be written as following: 

(X 1 e 1 \v(r)\X' 1 e[) J j = (jdMrMj^ , (2-36) 

and the particle-hole matrix elements: 

<Atiai|«(r)|Miai>^ = (hk\v(r)\j 3 j 4 )i r ■ (2-37) 

The matrix elements of Eq. (I2-37D which generate the interaction with the core 
clusters are not considered by the microscopic calculation of Refs. 6, 9) where the 
effective Hamiltonian is obtained by summing only over ladder diagrams. 4 ) By using 
the diagonal matrix elements of Eq. (|2-35p and the off-diagonal matrix elements 
given in Ref. 3) we diagonalize Eq. (|2-20p and obtain the amplitudes x( a i) an d 
x{ot2) of the mixed modes (2p) and (3plh). The calculation of the ground and 
excited state distributions and of the magnetic moments is then performed in terms 
of these amplitudes. The formula of the correlated distribution for the ground state 
of even nuclei is given in Appendix H. 



§3. Results 



In order to perform structure calculations, we have to define a single-particle 
base with the relative "single-particle energies" and to choose the nuclear two-body 
interactions. The single-particle energies of these levels are taken from the known 
experimental level spectra of the neighboring nuclei 27 ^ (see Appendix C) and given 
in Table I. For the experimentally unknown single particle energies of the fp shells 
we use the corresponding energies for the mass A-9 nuclei scaled accordingly the 
different binding energies. Our single particle energies agree reasonably well with 
those calculated in Ref. 28). Some of these levels are not bound. In this paper we 
perform as in Ref. 9) calculations by assuming all levels as bound. It has however 
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to be remember that the energies of Table I are used only in the first stage of the 
iteration procedure explained in Appendix F). 

For the particle-particle interaction, we use the G-matrix obtained from Yale 
potential. 29 - 1 These matrix elements are evaluated by applying the e s correlation 
operator, truncated at the second order term of the expansion, to the harmonic os- 
cillator base with size parameter b=1.76 fm. This value is consistent with the value 
used by Kuo. 30 ) As elucidate in Ref. 3) the potential used by the BDCM is sepa- 
rated in low and high momentum components. Therefore, the effective model matrix 
elements calculated within the present separation method and those calculated by 
are pretty similar. The separation method generates matrix elements which 
are almost independent from the radial shape of the different potentials generally 
used in structure calculations. 

The particle-hole matrix elements could be calculated from the particle-particle ma- 
trix elements via a re-coupling transformation. We prefer to use the phenomenologi- 
cal potential of Ref. 24). The same size parameter as for the particle-particle matrix 
elements has been used. 

One can generate the center-of-mass (CM) spurious states according to Refs. 32, 33) 
and evaluate the overlap between these states and the nuclear eigenstates of the 
model (see Appendix B). Model components having with the corresponding CM 
components an overlap greater than 10% were treated as spurious states and dis- 
carded. This is a convenient approximation considering that in our model space the 
energy of the CM is varying between 18 and 20 MeV. 

By using a base formed by 9 (2p) and 112 3plh-states we calculate the dis- 
tribution of the ground state of 6 He given in Fig. [TJ The spectroscopic factor for 
the ground state wave function, defined in Appendix G, and the most significant 
components of the ground state wave function are given in Table IX. 

In Fig. Q]we plot three distributions: 1) the correlated charge distribution of 6 He 
calculated with Eq. (|F-6p . 2) the correlated charge distribution of 6 He calculated with 
Eq. (|F-6P but neglecting the folded diagrams, 3) the charge distribution calculated 
for two correlated protons in the lsi shell. A charge radius of 2.25 fm has been 

2 

obtained for the distribution a), a radius of 2.39 fm for the distribution b), and a 
radius of 2.09 fm for the distribution c). 

In Table (II) the calculated charge radii are compared with the radii calculated 
by the other theoretical models and with the radius obtained by the IS theory. The 
radius calculated with the distribution c) (single particle) is in agreement with the 
radius obtained by the IS theory and with those of the other theoretical models. 
The BDCM increases however the charge radius of 6 He. One reason for this slightly 
disagreement can relay on the effect of the correlation operator which is neglected 
either in the evaluation of the charge radii via the Is theory or in the other theo- 
retical model calculations. As shown in Ref. 3) the correlations are also important 
in the analysis of the matter radius of 6 He. Here the larger matter radius repro- 
duces reasonable well the proton elastic scattering cross section measured at GSI at 
717 MeV/u calculated by using the Glauber method in the whole range of data. 35 ) 

The energy of the low lying spectrum of 6 Li is calculated by defining a coupled 
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base in which a proton-neutron pair in the (psd) single particle scheme of Table (I) 
are interacting with the (3plh) states generated by exciting the hole from the lowest 
(Is) shell. Within the model dimension of 13 (2p) and 196 (3plh) components we 
define the ground state wave functions. The spectroscopic factor for the ground 
state of 6 Li and the most significant components of the ground state wave function 
are given in Table X. By using this model space we obtain the spectrum shown in 
Fig. [2] right. For the ground state this model space is however not large enough to 
reproduce the second 1 + level and the magnetic moment of the ground state 36 ) . 
The magnetic moment of the ground state calculated in this model space is 0.875 
nm therefore larger then the experimental value (0.822 nm). The second spectrum is 
calculated introducing for the ground state a much larger space with dimension of 525 
components (13 two particle shell model states and 512 (3plh) states) which includes 
also the excitation of the particles from the \p-t state. Within this dimensional space 

2 

the spectroscopic factor for the ground state of 6 Li, defined in appendix G, and the 
most significant components are given in Table XI. Within this space we reproduce 
also the energy of the second 1 + level given in Fig. [2] left. To be considered that in 
other theoretical models the 1 + is always lying at too high energy. Since we plot 
the spectrum relative to the ground state energy, the effect of enlarging the base 
also for the other states is negligible at least for the first plotted 2 + and 3 + levels. 
The magnetic moment of the ground state calculated for the large model space is 
0.821 nm, value that is in good agreement with the experimental value. The charge 
distributions calculated within these two model spaces is given in Fig. [3j The two 
distribution are labeled: 1) charge distribution calculated for the small configuration 
space and 2) charge distribution calculated for the large space. The calculated charge 
radii are respectively (^h) 5 e< l ua l to 2.627 fm for the small model space and to 2.55 
fm for the large space. Both values reproduce well the charge radius of 2.55 fm 
obtained in Ref. 17) from the electron scattering experiments. Calculation of the 
electric scattering form factor measured in Ref. 17) for both charge distributions is 
presently under consideration. 

§4. Conclusions 

In this paper we have investigated the effect of the microscopic correlation op- 
erators on the spectra and the charge distributions of 6 He and 6 Li. The microscopic 
correlation has been separated in short- and long-range correlations according the 
definition of Shakin. 2 ) The short-range correlation has been used to define the effec- 
tive Hamiltonian of the model while the long-range is used to calculate the structures 
and the distributions of exotic nuclei. As given in the work of Shakin, only the two 
body short range correlation need to be considered in order to derive the effective 
Hamiltonian especially if the correlation is of very short range. For the long range 
correlation operator however the "three body component" of the correlation operator 
is important and should not be neglected. The effect of the "three body correlation 
operator" is to introduce in the theory a three body interaction. Therefore the use 
of the genuine three body interaction of the other theoretical model could, in the 
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present theory, generate double counting of diagrams. 

By using generalized linearization approximations and cluster factorization coeffi- 
cients we can perform expedite and exact calculations for the structure of 6 He and 
6 Li. Very good results for the spectrum of 6 Li have been obtained by considering 
large configuration spaces. Within correlated distributions we obtain charge radii 
slightly larger than those calculated within either non correlated distributions or IS 
experiments. 

This result should serve as motivation to a reevaluation of the different terms of the 
IS theory: i.e. the MS, which is mainly calculated in a perturbative approximation, 
and the FS, which is generally calculated in a point nucleus approximation, should 
be reevaluated within the non perturbative BDCM. 

Appendix A 

Definition of the model CMWFs 

In the BDCM the degree of linearization applied to the commutator equations 
defines the CMWFs of the model. For A=6 the model space is formed by two 
valence particle states and by the full set of the (3plh) CMWFs. These different 
components are associated to the following linearization mechanism: a) In the zero 
order linearization approximation we retain only two particle states: 

# 2p tiihJ) = la]A/ M \°) (A-l) 

For the two particles we distinguish between : 

1) effective valence space which is used to diagonalize the EoM, 

2) complementary high excited single particle states which are used to compute the 
G matrix. 

b) In the first order linearization approximation we include in the dynamic theory 
also the (3plh) terms. These are generated by the application of the correlation 
operator of the third order to the particles in the open shell states. Within this 
linearization approximation the CMWFs of the model are defined by: 

^resse d{jij2j) = + [a] 3 % J ) JM | 0) . (A-2) 

The (3plh) CMWFs are then expanded according to Eq. (|2-22|) . This expansion 
allows to orthogonalize the CMWFs in an easy way. c) The (4p2h) states which 
characterize the second order linearization step are not included in the model space 
but, linearized, generate the eigenvalue equation of the model (2p)+(3plh) states. 

Appendix B 

Center of mass correction 

Before performing the diagonalization of relative Hamilton's operator in the 
CMWFs defined in appendix A) we have to eliminate the spurious center of mass 
states. We start to compute, following the calculations of Refs. 32, 33), the percent 
weights of spurious states in the model wave functions. These can be obtained by 
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calculating the energy of the center of mass according to the following equation: 

e r = J dR^ dressed (jijj j) (i? 2 )if rdressed (iiij J) 

+2^| dndr 3 ^ dressed (j t j 3 J){n ■ f^^ti'^J). (B ' 1} 

»j 

In Eq. (|BTj) the calculation of the integrals can be performed by using the expansion 
for the (3plh) states given in Eq. f|2-22j) and by considering that for two particle states 
we have: 

(jdj J \{n ■ fj)\jijjJ) 

-*[**]( 4 i |) 2 {U ?}<w. (B - 2) 

where: 

j = (2j + 1). (B-3) 



By diagonalizing the above operator in the model space we obtain the energy of the 
center of mass. The overlap with the model space give the degree of "spuriosity" of 
the different components. 



Appendix C 

The single particle energies 

The model space which characterize the BDCM is formed by adding even particle 
to a closed-shell nucleus. The closed shell configuration can be described by a single 
Slater determinant and one can use the Hartree-Fock's theory to obtain the binding 
energy and the single-particle energies. Alternatively one can remark that for a 
closed shell nucleus (Z,N) the single particle energies for the states above the Fermi 
surface are related to the binding energies differences: 

e> = BE{Z,N) - BE*(Z + 1,N), (C-l) 

and 

e> = BE{Z,N) - BE*(Z,N + 1). (C-2) 
The single particle energies for the states below the Fermi surface are given by: 

e< = BE*{Z -1,N) - BE(Z,N), (C-3) 

and 

e< = BE*{Z,N -1) - BE(Z,N). (C-4) 

The BE are ground states binding energies which are taken as positive values, and 
e will be negative for bound states. (BE* = BE — E x ) is the ground state binding 
energy minus the excitation energy of the excited states associated with the single 
particle states. Within this method, which recently has been reintroduced by BA. 
Brown, 37 ) we derive the single particle energies from the known spectra of neighbor 
nuclei (see Table I). 
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Appendix D 

Iteration procedure to calculate dressed eigenstates 



Let us suppose we have two particles in the sd shell model states interacting via 
the (3plh) CMWFs. The configuration mixed wave functions are: 



i7>dressed 
w 2p 



l lrf 5 a 2si. 



+ [a 



t 

ld 3 



+ 



3 [4s i ais 1 } 1 

5 $ 1 



+ 



l L3 a 2sjVi rf5 aisi] 3 

2 2 7 3 



1 2 



(D-l) 



= 01 + 02 + 03 + 04 

solution of the eigenvalue matrix: 



(if) 2 \v(r)\ipi) 
{ip3\v(r)\ipi) 

V (04KOI0i) 
■ Vilo) \ 

V> 2 |o> 

Ho) 

4 |o) J 



(0i|u(r)|02) 
^2 + (^ 2 |v(r)|V 2 ) 

(0 3 Nr)|0 2 ) 
(^4^(0 1 ^2} 



(^i|u(r)|T/'3) 

(0 2 K r )l03) 
£3 + (03 KOI 0s) 

(04 KOI 03) 



(0iKOI04) \ 

(0 2 KOI04) 

(ip 3 \v(r)\ip4) 
£ 4 + (04KOI04) J 







(D-2) 

By diagonalizing the matrix of Eq. (jD-2|) we obtain four eigenvalues E\, E2, E3, E^ 
and the four eigenvectors given below: 

^1 = Xil0i) + X2I02) + Xsl^3> + xll04) (D-3) 
^2 = X?|0i) + xl\H + xl\H + xl\H (D-4) 
^3 = xflV'i) + xl|02) + + x!l04) (D-5) 

^ = xfl0i) + XalV-2) + x 3 |0 3 ) + X\\H (D-6) 

This eigenvalues and eigenvectors are then used to diagonalize the eigenvalue matrix 
in the second iteration step: 

/ E 1 +J* 1 \v(r)\¥ 1 ) (^Kr)|# 2 ) (^iKOI^s) (^i\v(r)\^) \ 

(^KOI^O £ 2 + (^ 2 K0l^2) (^3|«(r)j# 3 ) (^KOI&i) 

(l^ 3 |«(r)|^2> .Ba + ^slt'WI^) (^KOI&i) 

(^4|«(r)|^ 2 ) («i|«(r)|# 3 ) £ 4 + (^4K0l^4) / 

(D-7) 



V 



(•^KOI^i) 
(^KOI^O 



where: 



(*i|«(r)|!?i) 

= (Xi) 2 (0lKOI0l) + (X2) 2 (^2k(r)|V 2 ) + (Xs) 2 (03KOl03) + (xl) 2 (04KOl04) 
+2XiX 2 (0iKOI02) + 2x1X3 (01 K0 103) + 2xixi(0lKOI04) 

+2X2X3 (02 K r ) 1 03) + 2X2X4(02^(^)104) + 2X3X4 (03 KOI 04) 

(D-8) 

where the x}' s are the projections of the truncated model space on the basic vec- 
tors 2p, 3plh. The procedure is re-iterated until the energy convergence has been 
obtained. 
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Appendix E 

Basic equations for the CFT of the (Splh) CMWFs 



The normalization factors for the it operators are calculated by using the re- 
coupling algebra of Ref. 26). We obtain: 



coefi= y v[kj](-iy +t ' +k+j2+J+1 



J r k J\ 
J J 2 J r 



coef 2 = y (i;^)V[«/i«/ 2 «/ J i r l J s i^j](-i) i ' +j;+J ' l+J '' +J+1 

k 3 y 3 3 j* 3 i, 3 v 3 s 



(1 _§., 25jl{ _ 1) h + h-J i)§JsJ ^ 
(E-l) 



h 32 J\ 

J J 2 J'r 



h k h 

J'r jl Jr 



% j y j y 

J 32 Jr 



k J 



7 4 

s °r 
J J 



(i-«M-i(-i) il+i '" Ji ). 

and 

Coef 3 = V V[J l JprJskJ](-l) l+i ' + ^ +J r +J r 



X j y j ,y. 

i' k Js 

j j, J y. k 

(E-2) 



•<- E 

k J ^ 3^ 3^ 3 %3 y 3 s 



31 
J 



32 

J2 



Ji 
J', 



33 

J', 



3i 
31 



J2 
J} 



1 

J 



J y. 

31 



J y 



k 



Js 

J 



J y 

J 



I 

Ji 



j y 

k 

j y. 



(i-*Ms(-i) i!+jWi ). 

(E-3) 

By using these normalization coefficients we obtain the CFT for the three partitions 
by diagonalizing the following matrix: 



Js 

k 



Coefi*Coefi 


Coefi*Coef2 


Coefi*Coef3 


Coef2*Coefi 


Coef2*Coef2 


Coef 2 *Coef3 


Coef3*Coefi 


Coef3*Coef2 


Coef3*Coef3 



V 



33 n 



Js\J\ 



( b vf' 1 \a 2 J l J 2 \}X 1 e l Jre\J s )-Wf^(a 2 J 1 J 2 \}\ l e 1 J r e l J ! 



5^(3,1), 



J 

» j '{a 2 JiJ 2 \}\ 2 e 2 J r e 2 J s ) ■ 5 Vj 3 ' 1 \a 2 J 1 J 2 \}X 2 e 2 J r e 2 J s ) , . , , . ,. . ;.! , • , 

V 5 v} {3 ' 1 \a 2 J 1 J 2 \}X 3 e 3 J r el,J s )- 5 vi 3 ' 1 \ a2 J 1 J 2 \}X 3 e 3 Jre- 3 J s ) ) \ \([a)A j3 ] Jr K a i 

(E-4) 

The CFC coefficients associated to the u operators are derived within the same 
computational method introduced for the it operators; the normalization factors 
{cm ; fJ'i = 4, 5, 6) are given below: 

Coef 4 = Y v / [^i^i r 1 «/ s J J A;«/](-l) J+jl+J3+Jr2+J2+Jr+Js+J+1 



\a\a\ 



a' 3 a j4 



"11 



|J.)J) 
|J.)J) 
|J.)J) 



kj ' r 3^ 3^ 3 %•} V 3 ' s 



31 

J'r 



32 

k 



Jl 



33 

Jl 



34 

J 



J2 

J'r 



k 



Jl 



Jr 
J 2 



J'r 



32 



Js Jr 


J 


J 2 r 34 


k 


J'r 33 


J 


(E-5) 





Coef,5 



( Correlated EoM and. . . ) 
=Coef 4 * (-\)h+h--h with ( Jx 
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Coef 6 = yJ[kJ]{-l) ri+jA+l+k+Jl+Jr 
kJL, J f 



j 



k 
Ji 
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Ji 
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By using these normalization coefficients we obtain the CFT for the three partitions 
by diagonalizing the following matrix: 

/ 



Coef4*Coef4 


Coef 4 *Coef5 


Coef4*Coef6 


Coef5*Coef4 


Coef5*Coef5 


Coef5*Coef6 


Coef6*Coef4 


Coef6*Coef5 


Coef6*Coef6 



l(K«. 



/ 5 u\ {3 ' 1 \a 2 JiJ2\}^imJrViJs) - b uf l) {a 2 j 1 j 2 \}n im J r r}xJ s ) 
5 uf' 1 \a 2 JiJ2\}^V2JrmJs) ■ "'uf 1) ( a2 j 1 j 2 \}^mJrmJs) 



iil 
44) 



1 "■ Ur [a\al ]Js ^ 



J2 |3 



.<l';J Js ) J ) 



(E-8) 



Appendix F 

A numerical application of the CFT ■ 



In this appendix we apply the method of the previous appendix to calculate 
the CFC for a (Ids 2s i ) 2 (d:ipJ 1 ) 1 -(3plh) CMWFs formed by coupling two particles 

2 2 25 

assumed to be in the Ids and 2s 1 single particle shell model states to the dzp^ 1 p-h 

22 2 2 

pair. By using Eq. f|E- lj) we write for Coefi: 



Coefi 



Jr 



1 



(FT) 



which by restricting the quantum numbers to k=l for (Ji = 2, J 2 = J s = 1, J = 
l,J r = 2) give the cases of Table (III). The Coef2 given in Eq. (|E-2p assumes for 
k=l the form: 



Coef 



2 , £ 

J r JtJs 



(j;j r 4 )[5-3-9JiJiJ s ]2(-l) 



•j y 



3 
2 

'J ,. 



5 n 



2 J r 

1 3 



1 



1 




In Eq. (|F-2p the range of the indices in the sum is running over the following possi- 
bilities: 



J J> 


= 3/2 


J' r = 5/2 


J j> 


= 7/2 




= 0,1,2,3 


J}.= 1,2,3,4 


J l r = 


2,3,4,5 


Is 


= 0,1,2,3 


J s = 0,1,2,3 


Js = 


0,1,2,3 


J y> 


= 1,2,3,4 


J r 4 = 0,1,2,3,4,5 


J *. 


2,3,4,5 




= 0,1,2,3,4,5 


J r = 0,1,2,3,4,5 


J r = 0, 


1,2,3,4,5 


Ji 


= 1,2,3,4 


Ji = 1,2,3,4 


Ji = 


1,2,3,4 



(F-3) 
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By summing over all the possible cases we obtain for this special case the Coef2 
given in Table (IV). The Coef3 coefficients are given in Eq. (|E-3|) and for k=l we 
obtain: 



Coef 3 




I \ 2 w § | i w \ j$ jm r i 

1 1 J' r j\ J' r | Jl J\ J f J r j\ J r 

(F-4) 

For this special example the calculated Coef3 are given in Table (III). We see that 
the matrix we have to diagonalize in order to get the cluster coefficients is of the 
order of eight. By introducing the coefficient of Table (III, IV) in Eq. (IE-4j) and by 
diagonalizing the derived matrix we get the V-CFC given in Table (V). Analogous 
calculations can be performed for the u operators. The Coef4 given in Eq. (|E-5j) 
assumes for k=l the form: 

Coef 4 = J] [J 1 J 2 J}JsXU} 1 H-l) i+jl+j3+J?+j2+Jr+Js+J+1 

k J r t/i cJ? J %J r J s 




I | 2 I | I I 1 1 | 2 2 1 Ji 2 L J 2 °r <*r U j2 3 1 



J.9 1 
1 
1 



(F-5) 

By summing over all possible partitions we obtain the Coef4 coefficients which are 
given in Table (fH|) . For the Coefs we use the previous formula by replacing j\ — ¥ j2, 
and we calculated coefficients given in Table (VI). The Coefg are calculated form 
Eq. (IE-7P which in this special example takes the form: 

Coef 6 =^3.*(-l)^{ \ I I }. (F . 6) 

Now since J r = 1,2 we have for Coefe the coefficients given in Table (VII). By 
using the coefficients Coefj,i = 4,5,6 we derive the CFC for the U operators given 
in Table (VIII). By recalling that the pair coupled to J r is active we write for the 
matrix elements calculated in the ((ds si ) 2 (c?3p^" 1 ) 1 ) 1 CMWF the following value: 

2 2 2 j 



1 \ 1 
3 a 



(d5Sid3p 3 1 \v(r)\d5Sid3p 

2225 222 2 

= (0.0617) 2 (d5Si\v pp (r)\d5Si) 2 a + (0.0589) 2 (d3 si \v pp (r)\ds si ) 2 
+(-0.2088) 2 (d 3 si\v PP (r)\d 3 si)l + (0.2352) 2 (<L si |v P p(r)|(Lsi) 2 
+(0.6029) 2 (d5d3 \v pp (r)\dsd3) 2 a + (-0.3346) 2 (d5c4 \v pp {r)\d k d3 ) 2 a 
+ (0.6466) 2 {dld\\v pp (r)\dld\)l + (0.0908) 2 (s ips \v ph (r)\sip3_)\ 
+(-0.1504) 2 (sip3|^(r)| 2 S ip3)2 + (0.4022) 2 (d|k|t; ph (r)| 2 ( i3 2 p|)a 
+(-0.3278) 2 (dJpJ|^(r)|/|p 2 3) 2 + (0.7873) 2 (/|/| I^WIdJ/s) 2 
+(Q.2m^) 2 {dr^3\v vh (r)\d 1 p3)\ + (0.1883) 2 ((i5p3|L /l (r)|d5p 2 3) 2 ? 



(F-7) 
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One has to note that in the ladder approximations only the first seven terms of 
Eq. (|F-7p contribute to the matrix elements. 

Appendix G 

Spectroscopic factors of the dressed wave functions 

The spectroscopic factor of the ground states of mass A=6 isotopes is defined 

by: 

S P3P3 = (^14 3 flj, 3 |0) (G-l) 

The spectroscopic factor for the two neutrons in the 6 He together with the more 
significant components of the dressed J=0 + , T=l wave function of the ground state 
are given in Table IX. For 6 Li the spectroscopic factor and the more significant 
components of the J = 1 + T=0 ground state wave function are given in Table X 
for the small configuration space. Corresponding values calculated with the large 
configuration space are given in Table XI. 

Appendix H 

Charge distribution for two correlated particles 

In this appendix we give the difference between the distribution calculated for 
two non-correlated particles in shell model states and that calculated for two dressed 
(correlated) particles. In the shell model the distributions are evaluated by using 
the expectation values of the operators : 

P( r ) = '^2( a \ s ( r ~ r a )\a)ala a (H-l) 

a 

between the two particle shell model states: 

[(jii 2 ) J ]=a] i a J t 2 |0). (H-2) 

By performing small algebra we obtain: 

p{r) = p jl (r) +Pj 2 (r), (H-3) 

where pj i is the single particle distribution. This distribution is valid only in an 
extreme single particle model i.e.: Shell Model, Hartree-Fock, mean field theories. 
In the BDCM the effect of the long range correlation must be included consistently 
in the calculation of the densities. We need therefore to calculate the distribution 
starting from the correlated particle pair given below: 

#12 = e 5 ^i2 = (1 + Si + S 2 + S 3 + • • • )#i 2 , (H-4) 

where the Si,i = 1, 2 • • • are the correlation operators of the i-th order. Within our 
approximation the density should therefore be calculated from the model CMWFs: 



^12 = Yl xij^ij + Yl Xijki^ki- 

ij ijkl 
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By using the amplitudes of Eq. (12-21j) and the CFC of Eq. (12-221) we derive, discarding 
the isospin quantum numbers, the following equation: 

p J (R) = (J2xl 1 Y (^i^i^2^A|}n/iVLA) 2 angcoef(a 1 aiAAJ)^ i ( J R) 

01 nlNLX 

+2 Y &fc£ 5 ^ 3,1) («2JiJ2|}^ r eV s ) 

Y (ni/in 2 /2A|}n/A^LA)(n i /i</-A / |}n/iVLA / )angcoef(aie 4 AA / J r )^ L ( J R) 
nlNLXX' 

+ Y b uf 1 \a 2 JiJ2\}aiJra l J s ) 

Y (ni/in 2 /2A|}n/A^LA))(n,/ i n^iA / |}nWLA / )angcoef(aiaiAA / J r )^ L ( J R)) 

nlNLXX' 

+ Y Xa 2 X/3 2 ( Y 5 V} t(3 ' 1) (a 2 JiJ 2 |}eiJrei^)Vj 3 ' 1) (/3 2 JiJ 2 |}e i J;e, J s ) 

Q2/32 J r J r ieiej 

Y {n i lin[l' i \\}nlNL\)(n j l j n'/ j \'\}nlNL\^ 

nlNLXX' 

+ Y 5 uf 3 ' 1 \a 2 J 1 J 2 \}p i J r PM 5 U^ 1 \p 2 J 1 J 2 \}P j J' r P J J^ 

J r J r l0il3j 

Y (^/in-/-A|}n/iVLA)(n j / j n^^A / |}n/iVLA , )angcoef(/3^ j AA'J)^ L (^)))^ 

nlNLXX' 

(H-6) 

where R is the center of mass of the dressed particles, (nlNL) the relative and center of 
mass angular momenta of the ej) pairs , and the brackets are the Moshinski's 

brackets. 38 )> 39 ) The Xai are the shell model amplitudes and the Xa 2 the (3plh) 
amplitudes. The angcoef's of Eq. (|2 37|) are angular momentum transformation 
coefficients between the (jj) and the (Is) coupling given below: 



angcoef^aiAA'J) = hl 2 X\}^h)ji, {\h)h, A)(i± S ', W\}(¥x)&, (sM. A 
= V[JU2]A(-1) A <( h l 2 A 



i i S , , 

, 2 , 2 ; f L Z A 



Ji J2 J; 
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S J\2 J 



j\m\i(-i) x ' { a v 2 \' 

Jl J 2 Jl>2' 



a 3 f, f l r a' 



?' V T 



(H-7) 

The angcoef(aieiAA' J), angcoef^ie^AA' J), and angcoef(/5i/^AA'J) coefficients have 
a form analogous to that of Eq. (|H-7|) . The symbols a±, a'±, {3i,(3[, and ei,e[ are 
given below: 

«i — > two valence particles 

and (5[ — ► two particles from the (3plh) CMWFs (H-8) 
e\ and e[ — ► particle-hole form the (3plh) CMWFs 

The effect of the folded diagrams on the calculated distributions is given in Fig. Q] 
where we compare the correlated distribution of 6 He calculated with Eq. (|H-4p by 



( Correlated EoM and. . . ) 
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neglecting the particle-hole diagrams with that calculated by including also these 
diagrams. Analogous calculations will be performed for any operators. This effect 
has been until now not considered by the other theoretical models. The translational 
invariant density of Ref. 40) is derived by assuming that the wave functions of the 
nuclei are given in the terms of non correlated Slater's determinants. 
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Table I. Single-particle scheme and single particle energies (MeV) used to form the model CMWFs 



charge radius of 6 He Model 

1.944" 34 ' no-core shell model 

2.09 9 ' quantum Monte Carlo technique 

2.25 this work BDCM 

2.39 this work BDCM without the folded diagrams 

2.06 this work two correlated lsi —protons 

1.99 n) Cluster 

1.99 12) Cluster 

1.99 13) Cluster 
2.054 ± .014 10) Isotopic Shift (Exp.) 

Table II. Calculated charge radii for 6 He in fm compared with the results obtained in other theo- 
retical models and with the radius derived within the IS theory. 
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Table III. Coefi as functions of J r , Ji 
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Table IV. Coef 2 and Coef 3 as functions of J r , J. 
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Aj J r 



CFC 



0.0617 
0.0589 
-0.2088 
0.2352 
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-0.3346 
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Table V. Cluster factorization coefficients Vi ' ([ds si Ji = 2][d 3 p3 1 J 2 = i\}\jijjJr\\jij m 1 Js]) cal- 
culated for the seven allowed partitions and related CMWFs. 
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Table VI. Coef4 and Coefs as functions of J r , J. 
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Table VII. Coef6 as functions of J r , Ji 
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Table VIII. Cluster factorization coefficients U^' 1 ([d&s 1J1 = 2][d 3 p3 1 J2 = 1 Jr] [ji jm Js]) 

calculated for the seven allowed partitions and related CMWFs. 



Spect. fact. Ips lpa lpTlpT (Ip3ld5) 27r (lp3 ls^ 1 ) 2 '" (lp s Ids f' 1 (lps ls7 



tvtt 



0.9370 



0.9680 



0.1816 



-0.1145 



0.0628 



Table IX. List of the most significant components of the ground state + ,T=l wave function of 
6 He with E=-24.97 MeV 



Spect. fact. lpilpa lp 3 lpi lpilpi ld&ld3 (lps Ids ) 2 ' u (lp 3 lsi 



-1\1,U 



(lpi2si) D;CT (lp3lSi 



HTTTJ 



0.7649 



0.8746 



0.4559 



.1009 



0.0385 



0.0261 



-0.0111 



Table X. List of the most significant components of the ground state 1 + ,T=0 wave function of 6 Li 
calculated within the small base with E=-19.30 MeV 
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Spect. fact. TpHpl TpTTpl ld s Ids Ids Ids (lpa Ids )^ u (lp3 IsT 1 ) 1 '" (lpi2s i ) u ' u (lpi lsT 1 ) 1 '' 

2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 

0.6427 0.8017 0.4906 -0.0636 -0.1155 0.0653 0.03725 

Table XI. List of the most significant components of the ground state 1 + ,T=0 wave function of 6 Li 
calculated within the large base with E=-21.2 MeV 
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2',1 2M 2M 

2 ,0 2 , ,0 2 ,0 

0*,1 0*, i 1 

3*,C 3 ,C 3*,0 

v,o i*„q r,o 

Exp. BDCM BDCM 

with p-hole without p-hole 



Fig. 2. Calculated spectrum of Li: Left: energy levels calculated by allowing the excitation of the 
lsi- and the lp3-hole. Right: energy levels calculated by restricting the hole to the lsi. 
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radius (fm) 



Fig. 3. Calculated charge distributions of Li: black- calculated by restricting the hole to the si- 
state; white- calculated by considering si - and pj_-hole. The configuration space considered is 

2 2 

formed by 13 two-particle (shell model)-states and 512-(3plh) states. 



